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ABSTRACT 

Laboratory measurements of unpolarized and polarized absorption spectra of 
various samples and crystal stuctures of silicon carbide (SiC) are presented from 
1200-35,000 cm"^ (A ~ 8-0.28 /im) and used to improve the accuracy of optical 
functions (n and k) from the infrared (IR) to the ultraviolet (UV). Comparison 
with previous A ~ 6-20 /xm thin-film spectra constrains the thickness of the films 
and verifies that recent IR reflectivity data provide correct values for k in the 
IR region. We extract n and k needed for radiative transfer models using a new 
"difference method" , which utilizes transmission spectra measured from two SiC 
single-crystals with different thicknesses. This method is ideal for near-IR to 
visible regions where absorbance and reflectance are low and can be applied to 
any material. Comparing our results with previous UV measurements of SiC, 
we distinguish between chemical and structural effects at high frequency. We 
find that for all spectral regions, 3C (/9-SiC) and the El.c polarization of 611 (a 



^currently at Planetary Science Institute, Tucson, AZ, 85719 
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type of a-SiC) have almost identical optical functions that can be substituted for 
each other in modeling astronomical environments. Optical functions for i?||c of 
6H SiC have peaks shifted to lower frequency, permitting identification of this 
structure below A ~ 4/im. The onset of strong UV absorption for pure SiC occurs 
near 0.2 /xm, but the presence of impurities redshifts the rise to 0.33 /zm. Optical 
functions are similarly impacted. Such large differences in spectral characteristics 
due to structural and chemical effects should be observable and provide a means 
to distinguish chemical variation of SiC dust in space. 

Subject headings: Silicon carbide - Optical constants - Methods: laboratory - 
Stars: carbon - (ISM:) dust, extinction - Infrared: stars 



1. Introduction 



Silicon carbide (SiC) is of great interest in the context of dust in a variety of astro 
physical environments and is t he most well-studied type of presolar grain (IClayton fc Nittler 
2004J : iBernatowicz et al.ll2006l . and references therein). Isotope measurements exist for thou- 
sands of presolar SiC grains (mostly of the cubic (3 structural form) and suggest that they 
came from AGB stars and supernovae; the majority of meteoritic SiC grains are thought 
to have originated in circumstellar shells of carbon stars (C-stars). Astronomical SiC was 
discovered via observations of a broad ~ ll/xm emission feature seen in many C-star spectra 



(jHackwelllll972l : iTreffers fc Cohen! Il974j ) that is associated with stretching of the Si-C bond. 
In the absence of strong infrared (IR) features from amorphous or graphitic carbon grains, 
the presence and variations of the 11 /im SiC feature provide a useful diag nostic for inferring 



the physical conditions in the circumstellar dust shel l s of C-stars (e.g., iBaron et al 
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tion, evolution and destruction in many astronomical environments. [Thompson et al.l (120061 ) 
showed that relating observational data on the SiC peak position, strength, and shape to 
the underlying continuum temperature did not yield useful information about circumstellar 
shell structure and evolution, and concluded that radiative transfer modeling is necessary to 
disentangle the dust shell parameters. The ability to perform and achieve the latter has been 
limited by a lack of grain-size-independent optical functions (real and imaginary indices of 
refraction n(A) and k{X)) for SiC at all relevant wavelengths. 

Studying SiC presents several opportunities for understanding cosmic dust. In particu- 
lar, the prevalence of SiC grains around carbon stars and in presolar samples suggests that 
some SiC grains survive the interstellar medium (ISM) and are incorporated into new plane- 
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tary systems (e.g.. lBernatowicz et al.ll2006l . and references therein). Ho wever, the abuiidance 



of SiC in the ISM determined in most studies is extremely low (e.g., iWhittet et al.lll990 
Kemper et al.l 120041 ). with only one study e stimating abundances consistent with 10% of Si 
atoms being bound in SiC (IMin et al.l 120071 ). Because UV extinction is the main constraint 
in m ost ISM dust grain models, this disparity is connected to the high photon fl ux in the 
UV feim. Martin fc Hendrvlll994l : IZubko et all 120041 : IWeingartner fc DrainelboOlh . coupled 
with the fact that optical strengths f or SiC in the UV vary among published studies (cf. 
PhihpplEgSsI : Ichovke fc Patrick! 1 1 95 7h . 



In addition, recent Spitzer Space Telescope observations have shown that carbon stars 
are much more common in lower metallicity Local Group galaxi es than in our own Galaxy 
dlagadec etaPboO?! . boosi : Iziilstra et aPbood : Isioan et al.lboosh . Detection of an SiC fea- 
ture near 11 /xm provides a diagnostic for the nature of these objects. Furthermore, SiC is 
commonly fou nd in extreme carbon stars and planetary nebulae (PNe) in the Large Mag- 



ellanic Cloud (IStanghellini et al.l 120071 : iBernard-Salas et al 



not commonly found in their Milky W ay counterparts (IBarlow 
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Gruendl et aPboOsh and 



1983 



Rinehart et al. 2002 



Casassus et al.ll200ll : ISpeck et al.ll2009l ). Consequently SiC presents a tool for understanding 
the effect of metallicity on dust formation, evolution, and survival. Thorough investigation 
of these environments requires reliable optical functions for SiC with a broad wavelength 
coverage. 



Pitman et al.l (12008 



hereafter designated Part I) provided new mid- and far-IR re- 
flectance spectra and grain-size-independent optical functions for several SiC samples and 
discussed the limitations of existing mid- to far-IR datasets on n and k for SiC (e.g., use 
of the commercially available hexagonal a rather than the space-relevant cubic (3 structural 
forms of SiC, dependence on grain-size distribution of original laboratory sample measured, 
and the relationship between structure and spectral parameters from sy mmetry analysis). In 



particular. Part I conflrm ed that previous SiC reflectance studies (e.g.. ISpitzer et al.lll959 



Bohren fc Huffmanlll983l ) correctly assessed peak strengths and widths but found that peak 
positions were about 4 cm~^ too low, due to the low resolution and lack of internal cali- 
bration in the older instruments used. Quantitatively analyzed IR reflectivity data (Part I) 
and previous laboratory studies of SiC, particularly at high frequency, provide much of the 
information required for the applications discussed above, but improvements are needed in 
the following areas: 



1. Pubhshed visible-UV spectroscopic data for /5-SiC cover only part of the needed spec- 
tral range, and past visible-UV reflectance and absorbance datasets for a-SiC are not 
entirely consist ent with one another. For both a- and /3-SiC, n was obtained only 
for the visible JThibaultl Egii : IPhihppI Il958l : IPhilipp fc Taftlll96ol : IChovke fc Patrick 
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19571 ). Furthermore, those data are not in electronic form and have been incorporated 
into previous derivations of optical functions via digitization, which provides limited 
spectral resolution and accuracy. 

2. For SiC, the near-IR to visible region has low reflectivity {R < 20%) and low absorp- 
tion coefficients {A < 0.1 yum^^). As a consequence, not only do primary reflections 
occur at the front surface, but moreover, secondary reflections at the back surface of 
the sample are significant (Fig. [1]). These back-reflections affect both reflection and 
absorption/transmission measurements, thereby contributing errors to the extracted 
optical functions for these regions. 

3. In the UV region, both R and A are high, from which accurate optical functions can 
be extracted; however, few modern laboratories are equipped for measurements above 
about 33,000 cm~^. Very high frequency data in early literature can be evaluated and 
utilized if data are corrected for reflectivity are gathered in the visible. 



To address these issues, this paper provides new, single-crystal absorption spectra of 
the same samples used in Part I, but at higher frequencies and over the spectral range where 
R is low, to better constrain k up to the UV. Combining t hese new data with the o ptical 
functions from Part I and our previous thin-fllm IR data (ISpeck et al.l Il999l . 120051 ) both 
augments existing data by signiflcantly extending the wavelength coverage, and provides a 
cross-check on the accuracy of the previous data. 

In addition we develop a new approach - the "difference method" - for extracting n and 
k val ues from tran smission spectra and also utilize "limiting values" as in previous literature 
(e.g. iLipsonI Il960l ). Our approach can be applied to any material, not just SiC. Because 



transmission is easily and commonly measured, the difference method and limiting values 
may be used to greatly expand the database of astronomically relevant near-IR to UV optical 
functions and aid implementation of radiative transfer models for many environments, not 
just C-stars. 



2. Experimental methods 

We investigated seven different samples, listed in Table [H Although the prevalent form 
of astronomical SiC is the cubic (3C, or /3-SiC) polymorph, large single- crystals of that type 
are not available. Hence, we examined vapor deposited 3C as well as several crystals with one 
speciflc hexagonal, layered structure (6H). The label a is used to describe diverse hexagonal 
variations of the SiC structure. Part I describes sample characteristics in detail, which are 
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summarized in Table [H The colors of the samples are used to denote which sample is being 
examined. Because color is indicative of impurities, this attribute is significant, as will be 
discussed in § I4.3.1I 

The 6H SiC samples are anisotropic, so it is necessary to gather spectral data parallel and 
perpendicular to the layers (the a and c crystal axes) to provide data for the extraordinary 
and ordinary rays, respectively. Sample thicknesses were determined by using a calibrated 
micrometer and also by using a binocular microscope with a calibrated reticule. If possible, 
absorbance spectra were obtained from commercially polished sections. Thin samples were 
prepared by grinding and polishing sections of moissanite, a gemstone grade variety of 6H 
SiC. Twin planes provided reference points to maintain correct orientations (a-a or a-c 
plates). For the near-IR to visible studies of tiny platelets, thicknesses were obtained from 
interference fringes using: 



2nAz/ 



where d is thickness and Au is the wavenumber spacing of the interference fringes. Fringes 
tend to occur where d is similar to the wavelength of the incident light and when the faces 
are perfectly parallel. Values of n from our reflectivity study (Part I) and interference fringe 
spacing in the near-IR were used. Nearly identical v alues of d for thiii platelets were obta ined 



using direct measurements of n from refractometry (jShaffer &: Naunj|l969l : IShafferlll97ll ) and 
fringe spacing in the visible. 

Room temperature (18-19°C) IR absorption spectra were acquired at normal incidence 
using an evacuated Bomem DA 3.02 Fourier transform spectrometei0 (FTIR) at 1 cm~^ 
resolution below u = 4000 cm~^ and at 4cm~^ resolution at higher frequencies. Resolutions 
of this level are considered reasonable for non-gaseous samples acquired at room temperature. 
The accuracy of the instrument is ~0.01cm~^. Far-IR data {u < 650cm~^) were collected 
using a SiC globar source, a liquid helium cooled bolometer, and a coated mylar beam- 
splitter. Mid-IR data = 450-4000 cm~^) were collected using a SiC globar source, a 
liquid nitrogen cooled HgCdTe detector, and a KBr beam-splitter. Near-IR data (z/ = 1800- 
7000 cm^^) were collected using a SiC globar, a liquid nitrogen cooled InSb detector, and a 
CaF2 beam-splitter. For z/ = 4000-10,000 cm~^, we used a quartz lamp, a liquid nitrogen 
cooled InSb detector, and a quartz beam-splitter. To access the z/ = 9500-21,000 cm~^ range, 
we substituted a silicon avalanche detector. Large samples were mounted on apertures at 
the focal point in the sample compartment, and a wire-grid polarizer was used to obtain 
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polarization information from the moissanite sample. Most of our 6H samples were a-a 
plates, from which unpolarized data (equivalent to the orientation El.c, where E is the 
electric field vector) were gathered. Small platelets from Alfa/Aesar"""^ were placed on 
apertures of 100 yum diameter and data were collected using a Spectra- Tech IR microscop^. 

Additional, independent measurements were obtained from the platelets (on apertures) 
at UV wavelengths using the spectroscopic facihties at the Geophysical Laboratory, Carnegie 
Institution of Washington. This system consists of a custom-made, all-reflective microscope 
with Cassegrain-type mirror objectives and two off-axis paraboloids as a transfer optics, 
coupled to a Nicolet 750 Magna FT-IR spectrometer for IR measurements and to a grating 
monoch romator-spectrograph w ith a CCD detector, and configured for the visible and UV 



regions (iGoncharov et al.ll2006l ) 



We used a CaF2 beam-splitter and a halogen lamp for near-IR acquisitions, and a 
combination halogen/deuterium light source for the visible-UV spectral range. Spectral 
resolution was about 16 cm~^ through the whole spectral range (2400-40,000 cm~^). 



3. Data Analysis and the Difference Method 

Because transmission through a sample is more easily and commonly measured than 
reflection, we present two methods for extracting n and k from transmission measurements: 
the "difference" and limiting methods. Which method is applicable depends not only the 
material property A (or equivalently k) but also on the thickness. Because optical functions 
are frequency dependent, for any given thickness, each method will give accurate results 
over only certain spectral segments. Practical issues such as sample characteristics (Table [T]) 
strongly influence thickness and therefore which method is viable for the region of interest. 
In addition, the experimenter can vary thickness to optimize the accuracy of the results. 
The equations underlying both methods are definitions for the various optical properties. 
The following analysis assumes normal incidence of a coUimated light beam on a sample and 
that interference fringes do not exist or can be removed, e.g., by smoothing. Instruments 
record the light that is transmitted through a sample, rather than the light that is actually 
absorbed by it. Software for both spectrometers provide absorbance using the convention 
from the spectroscopy, chemistry and mineralogy literature (SCM), 

aSCM = - logio(^— ), (1) 
-'0 
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where Imeas is the intensity exiting the sample and /q is the intensity entering the sample 
(which is obtained by collecting a reference spectrum). The effect of surface reflections is 
neglected in the SCM literature largely because R is unknown. Instead, empirical baseline 
corrections are used, which suffice to extract widths and relative intensities needed for many 
applications. However, baseline corrections do not provide the information on R and n 
needed for modeling radiative transfer. 

To compute the amount of light that is actually absorbed, the true or ideal absorption 
coefficient (A) requires correcting ascM for specular (surface) reflections of magnitude R, 
where R is the proportion of incident light reflected only from the front surface as shown in 
Fig. [H as follows: 



Conversion to natural logarithms via Eq.[2]is needed to directly compare measured absorption 
spectra with the ideal absorption coefficients A obtained from reflectivity data (see, e.g.. Part 
I, Figures 3a-10a). 

Transmission spectra measured from two samples with different thicknesses, d, but that 
are otherwise identical, are described by Eq. [2] with different values for d. Taking the differ- 
ence between those equations (not shown) eliminates the reflection term, thereby allowing 
determination of A: 



A = 2.3026 ''*^'^' y^" . (3) 

dthck ~ dthin 

where the subscripts "thin" and "thck" indicate the two different sections being measured. 
Specular reflectance R is then obtained from: 



— I _ lQ{dthinA/2.3026-athi„)/2 



(4) 



We refer to this approach as "the difference method." To the best of our knowledge, this 
equation has not previously been published, although the above analysis is an extension of 
existing analysis in t he spectrosco py literature (i.e., use of limiting values for transparent 
spectral regions, e.g.. iLipsonlllQGOl ). described below. 



The imaginary part of the complex index of refraction, k, is defined as 

A 



k 



(5) 



8 



Since we can extract both k and R from the two transmission spectra, we can then obtain 
n, the real part of the complex index of refraction, from the following equation: 



which assumes that the surrounding medium (air or vacuum for single-crystal studies) 
has n = 1. 

Applying this analysis requires meeting certain conditions. The thicker section cannot 
be opaque over any part of frequency range being investigated. The thinner section cannot be 
so thin that its absorbance is negligible. The polish must be nearly perfect for both sections 
and the thicknesses must be such that fraction of light reflected from the front surface is 
similar to the fraction of light transmitted through each crystal. Also, any experimental 
uncertainties connected with measuring thicknesses will be propagated into errors in the 
optical functions. Multiple measurements can decrease the uncertainties associated with 
computing differences, provided that all sections are identical in composition and surface 
pohsh. 

Due to the stringent requirements for the difference method, we utilize cross-checks. 
One independent confirmation exists for n in the visible because in this spectral region, 
k -Cl for small to moderate values of A. For this case, Eq. [H] simplifies to 



which can be evaluated using measurements of n in the visible, obtained independently by 
matching the sample's optics to that of a liquid with known refractive index. 

A second cross-check, or alternative approach to our "difference method" is use of limit- 
ing values for determining n and k. For regions where A is nearly zero, spectral measurement 
of a fairly thin sample provides R because transmission, T = Imeas/h, is given by 




(6) 




(7) 



T 





1 - i?2e-2Ad 



and can be simplified for the case of negligible absorption to 



(9) 



- 9 - 



Substituting Eq. [7]into Eq. [9] yields n (e.g. lLipsonlll960l ). Both reflections from the front and 



back surfaces of the sample section are accounted for (Fig. [I]) . Use of limiting values requires 
a well-polished sample. To obtain A or fc, additional transmission data is then acquired from 
a very thick sample, and the results for R are used in Eq. [21 

If both sample thickness and A are very large, a special case occurs. When scaling 
large A obtained from 1 mm sections to thicknesses near 1 /im, the correction for reflectivity 
becomes negligible and A can be directly computed. Again, a good polish is needed or 
scattering occurs which has a frequency dependence. 

By applying both the "difference" and limiting methods, we can compute n and k from 
multiple transmission spectra of identical samples of different thicknesses. By comparing 
the results over frequency regions where both methods are applicable, we can construct a 
consistent set of optical functions and ensure that surface polish is not affecting the val- 
ues. Such a comparison requires laboratory measurements over a wide range of thicknesses. 
Furthermore, this approach can be applied to any material, not just SiC. 



4. Results and Comparison to Previous Studies 

4.1. Near-IR absorption of 3C (3-SiC 

Laboratory absorption spectra from a thin wafer of cubic SiC show one very strong 
IR mode from v ~ 950-1000 cm~^ that is off-scale (Fig. [2]). The noisy profile close to the 
barycenter is affected by stray light and is not intrinsic. This main feature results from 
Si-C stretching and is the transverse optic (TO) component. This triply degenerate mode 
also has a longitudinal optic (LO) component. Ideall y, no other IR mode is active in 3C. 



However, as is common in diatomic substances (e.g., iHofmeister et al.l 120031 ). structure is 
superimposed on this main band. A weak shoulder at 614.4 cm~^ is associated with the 
longitudinal acoustic (LA) mode, which is present due to resonance with the TO and LO 
modes. Transverse acoustic (TA) modes have frequencies too low for such resonances to 
occur. Peak positions of the fundamental modes are given in the footnotes in Table [21 

From u ~ 1000-1800 cm~^, a series of fairly intense overtone-combination bands exist 
that are well resolved for sample thicknesses of c? ~ 5-100 /xm. From v ~ 1800-2700 cm~^, 
weak overtone-combination bands exist that are well resolved for d ~ 100 /im to several 
mm. Positions of the overtones and the relationships to the fundamental modes are given in 
Table [H Most of the modes are clearly connected to some combination of the fundamentals, 
but some cannot be assigned with confidence, as indicated by the question marks in this 
table. It is difficult to predict the activity and frequency of overtone-combination modes 
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because thi s type of phonon involve anharmonic interactions of fundamental modes (e.g. 



Mitral Il969l ). Sometimes Raman modes can sum to provide overtones in the IR whereas 
some combinations of IR modes are forbidden. Overtones lacking an obvious assignment 
could be connected with impurities, non-stoichiometry, or stacking disorder. 

By comparing intensities of the shoulder near 1800 cm~^ from the spectra of the thin 
wafer to spectra of a much thicker wafer, we confirmed that the thin wafer indeed has a 
thickness of 25 fim (Fig. [2]). 

Absorbance spectra at liquid N2 temperatures of the 3C wafer were also acquired by 
attaching the thick wafer to a cold finger. No shifts in peak positions were seen. Bands are 
sharper but have roughly the same area as measured at room temperature (Fig. [2]). From 
this, we infer that absorption spectra of SiC are litt le affected by temperature, consistent 



with the incompressibility of SiC (e.g.. lKnittlelll995l ). i.e., bond lengths for silicon carbide 



vary little with pressure or temperature. Positions and widths of the main bands will change 
even less than indicated by Fig. [2] upon cooling because u of the fundamental mode changes 
at about half the rate of u for the simple overtones. 

The polish on the thin and thick sections were quite different (the vapor deposit has 
some texture) and opaque spectral regions exist, making it difficult to apply the difference 
method. For example, near 2500 cm~^, using Eqs. [3] and H] gives A as 0.00279 fim~^ and 
R as 0.048. This computation of R differs by a factor of 4 from the value in Part I and is 
attributed to variations in polish on the surfaces of the samples. We therefore determine the 
absorbance of the 1 fim wafers (§ I4.1.ip by scaling from the very thick (308 /im), strongly 
absorbing section. 



4-. 1.1. Composite spectra of 3C SiC; inference of thin-film thicknesses and band strengths 



Thin film spectra of bulk /3-SiC may not be precisely equivalent to spectra of the wafer 
because the bulk sample is not stoichiometric (ISpeck fc Hofmeisterll2004l ) and because it was 
difficult to create uncracked thin films of uniform thickness for this extremely hard material. 

We instead compare absor ption spectra of o ur wafer to thin-films of nano crystalline 
/3-SiC (Fig. [31 data shown from lSpeck et al.l 120051 ). In the thin-films, LO modes and shoul- 



ders are exag gerated due to a small proportion of the light crossing the film at non-normal 
incidence (see lBerremanlll965l ) . Because the refiections mainly arise from the smooth, highly 
polished diamond faces of the sample holder, rather than from the relatively rough thin-film, 
and are included in the reference spectrum, accounting for SiC refiections does not provide 
results that are consistent with refiectivity measurements and analyses (see results for MgO 
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by iHofmeister et al.ll2003l ). Therefore, we compared the raw thin-film data directly to the 
single-crystal measurements. We did not directly c ompare our data to available absorption 
spectra from dispersions (e.g.. lMutschke et al.lll999l ). because peak p rofiles are affected to an 



unknown degree by factors such as light leaking around grains (e.g. 
200l[ l. 



Hofmeister et al. 2000 



The thin-film data could be matched to the ideal absorption coefficient at either high 
or low frequency, but not at both. Fig. [3] therefore presents two possible merges for the SiC 
wafer and thin-film absorption spectra, one at 750 cm~^ and one at 1000 cm~^. Although 



inferred thicknesses for both me rges agree with results for other minerals (IHofmeister et al. 



20031 : IHofmeister fc Boweyll2006l ). the high frequency merge (suggesting an original thickness 
of 0.34 nm) is preferred because it gives more reasonable strengths for the overtones and 
is provided in electronic form (Table [3]). Overtones were not seen in thin-film spectra and 
must be sufficiently weak to be hidden by noise or interference fringes. A thinner wafer is 
needed to reduce the ambiguity in merging spectra, thereby more accurately determining 
the thickness of the thin film, but it is difficult to grind samples to the ca. 5 fim in thickness 
needed while maintaining parallel faces. Note also that our comparison of thin-film spectra 
to those of the wafer ignored the contribution of reflectivity as a function of frequency (see 

As a first approximation, the electronic file for the true absorption coefficient includes 
subtraction of 0.00145 /im~^ which was computed from the reflectivity correction term in 
Eq. 2 using R = 0.2 at 2500 cm~^. We neglect the frequency dependence of R because for 
the thick, strongly absorbing samples the correction term is small. 



4.1.2. Laboratory measurements vs. ideal absorption coefficients 

Our merged /3-SiC absorbance spectra can be reconciled with the ideal absorption spec- 
tra derived from classical dispersion analysis of reflectivity (Part I), as follows. Due to light 
leaking through cracks in the thin-film, the main peak and u ~ 800-950 cm~^ is rounded. 
The merged absorbance spectra "fits" into the ideal absorption peak at about half-maximum. 
Excess intensity near the LO position is expected, due to non-normal beam incidence and 
possibly a wedge-shaped film (the diamond faces in the sample holder are not perfectly par- 
allel). The excess intensity from the LO component also shifts the absorption maximum 
from the TO component towards the LO position, as seen in the comparison to the high fre- 
quency merged spectrum (Fig. [3]). This comparison supports the lower value of the thickness 
(0.34 /im), inferred through comparison in the previous section. 
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4.2. Near-IR absorptions of 6H samples 

Figure m compares measured absorbance spectra from 6H platelets of various thicknesses 
to a thin-film made from the same material. The thin-film may contain some particles 
oriented with -E||c, but this situation is unlikely because crushing would flatten the platelets 
rather than reorient them. A small (below 5%) component of ii^||c would not be detected. 

Two platelets with measured thicknesses of 5 fim gave slightly different absorbance 
spectra. Most of the differences are observed for the low intensity peaks and result from 
superposition of interference fringes on the weak peaks. Because the positions (spacing) 
and intensities of these interference fringes vary among the spectra taken from the three 
platelets, the weak overtone-combination modes can be discerned (Fig. H]). In addition, 
the far-IR data show some low frequency modes associated with zone folding (see Part 1; 



Nakashima fc Harimal 119971 ) 



Variations seen in the height of the main peak originate in slight differences in thickness 
that are within the uncertainty of the measurements. The fiat peaktops occur because of 
diffraction around the particles and because samples of this thickness are opaque at peak 
center, due to the combination of high reflectance and absorbance. Also, the peak height 
appears lower for the thick plate because of the combined effects of surface reflectance and 
stray light. 

We also thinned the sample of moissanite and obtained absorption spectra for both 
polarizations for several thicknesses (Fig. [5]). For the SiC platelets and moissanite, peak 
positions appear to be the same within uncertainties of 1-2 cm~^, depending on peak width. 
We therefore report peak positions in Table [2] obtained from the best resolved peaks in 
the moissanite spectra. Assignments are based on proximity of the peak positions obtained 
from summing the frequencies of the fundamental modes and on the shifts of peak positions 
between the two polarizations. 

Based on binocular microscope and micrometer measurements, the thickness of the 
thinnest moissanite piece is constrained to be 85-90 fim. In addition, the thickness of the 
thick platelet was uncertain because this sample was too small to measure using a microm- 
eter, and the determination using the binocular microscope was suspect. By comparing the 
average of the two absorbance spectra obtained from two thin SiC platelets to the absorbance 
spectrum obtained from the thick SiC platelet (Fig. H]), we infer that the latter had a thick- 
ness of 22 fim. Comparison with moissanite spectra gave a similar value. Note that Fig. H] 
shows a scaled spectrum, not the absorbance for a 22 fim section. 
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4-2.1. Cross-checks: comparing 6H and 3C SiC 

Hexagonal forms of SiC frequently have stacking faults and twinning that could produce 
extra, internal reflections and affect extracted values of R and A. We therefore compared 
measured absorbance spectra of 3C SiC (Fig. [2]) to those of our various 6H samples (Fig. H] 
and Fig. E]). We found that peak positions and spectral profiles differed only slightly among 
the various orientations and layering arrangements, and thus internal reflections are not a 
concern for the thinner samples. 

Similarity of the 3C to 6H spectra supports equivalence of spectra from the platelets 
and moissanite. The main difference in spectra occurs between the two polarizations of 6H 
and 3C largely results from the 10cm~^ difference between the TO peak positions of Si-C 
fundamental stretch in the two orientations. Intensities and widths seem to not be strongly 
affected by structural form. 

4-2.2. Composite spectra of 6H SiC; inferred thin film thicknesses and band strengths 

Spectral segments with the highest resolution and least noisy peaks were merged to 
form the composite absorbance spectrum for i?±c over all wavelengths (Fig. [6]), and for E\\c 
at high frequency. The merge suggests an original thickness of 0.54 /xm for the thin film 
and a thickness of 15 fim for the thinnest platelets. The former value is consistent with 
the thickness inferred for the thin-film of the cubic SiC sample, whereas the latter value 
is outside the uncertainty of the determinations made with the binocular microscope. The 
overestimation in the latter value largely results from the data containing specular reflections 
that effectively increase the absorption. The fit within the peak envelope of the classical 
dispersion analysis is good. If instead we assume that the thicknesses of the platelets are 
correct (i.e., we use an average of two thicknesses that were measured as near 5 /xm), then 
the original film was 0.18 /im thick and fits within the envelope except near the shoulders. 
The different thicknesses obtained for the thin film through comparison is largely connected 
with reflectivity not being subtracted. 

The merged spectrum for the 0.54 /xm-thick sample is provided in electronic form (Ta- 
ble Hj) for 6H because it is similar in thickness to the 3C film, and because we have only 
succeeded in making films as thin as 0.18 /im for soft materials, which SiC is not. For this 
spectrum, we spliced in segments in the far-IR from Paper I's classical dispersion analysis 
because our far-IR laboratory absorption data had strong fringes. Fringes are unavoidable 
because the sample thickness needed to resolve the far-IR peaks is optimal for creating fringes 
in this range (A= 25-100 /xm). The fringes are very similar to peak widths in this spectral 



-14- 



region which makes spectral subtraction or polynomial fitting difficult to perform accurately. 
Note that the fits are very close on the edges of the peaks. This spectrum can be multiplied 
by a factor of three, and still be within the uncertainties of our merge. However, the overtone 
region is best represented by the data as provided. For ideal values of A, the results from 
classical dispersion (electronic Table [5]) were substituted across the main peak. 

For the £'||c polarization (shown as light lines in Fig. [6]), merged data are provided in 
electronic Table O For the main peak, synthetic spectra from classical dispersion analysis 
(Part I) were used to produce an ideal spectrum. 



4.3. SiC absorptions in the visible to UV 

We studied a relatively thick 6H platelet using both the Bomem (FTIR) spectrome- 
ter and Nicolet-CCD spectrometer system. Data obtained from these two different instru- 
ments agree well for frequencies below 16,000 cm~^ (not shown). Above this frequency, high 
throughput an d high sensitivity of availab le detectors overcome advantages of using an FT 



approach (e.g., [Griffiths &: deHasethlll986l ) and so we focus on results from the CCD spec- 
trometer. 

Raw data obtained from two 6H platelets with surfaces as grown (Fig. [7]) contain in- 
terference fringes, the separation of which provides accurate thicknesses of 5.9 and 23.8 /zm, 
respectively. The absorbance of the thicker platelet is sufficiently high such that data ac- 
quisition was limited to u < 30,500cm~^. To avoid fringes in our derived optical functions, 
we fitted the results to 5th order polynomials above u = 20,000cm~^ (not shown). To 
extrapolate data for the thick platelet up to u =35,000 cm^^, we linearly fit absorbance of 
the thick platelet from 30,000 to 30,600 cm^^. This linear extrapolation provides the lowest 
possible values of A at the highest frequency attained, and the largest R values (Fig. [7 1) . The 



variation is small, however, and all such fits agree with A obtained by lPhilipp &: TaftI (Il960l ) 
through Kramers- Kronig analysis of their UV reflectance data. Data on the platelets were 
also fit to higher (up to 9th) order polynomials and extended to lower frequencies, but the 
results (Fig. [8]) were essentially the same. In all cases, calculated reflectivity is consistent 
with direct measurements and index of refraction measurements at the red end of the visi- 
ble, with a dip occuring near 29,000 cm~^, and a steep rise to high frequency. There is some 
leeway in the fitting, which apparently underestimates R and n values and overestimates A 
and k in the calculations near 30,000 cm~^, and vice versa elsewhere. 
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4-.3.1. Stoichiometry, impurity, and color effects 



In this subsection, we use the observed spectral differences and the known connection of 
impurities with absorptions in the visible-UV region to ascertain trends. Colors of samples 
generally indicated the presence of impurities. The colors of our samples are listed in Table [TJ 

The steep rise in A and R for the SiC platelets into the UV above p = 26, 000 cm~^ is 
due to the band gap in this semiconducting material. The spectral characteristics appear to 
be affected by non-stoichiometry and other types of impurities existing in t hese ma t erials , 
as indicated by coloration. For example, absorption m easurements of 6H by iPhilippI (119581 ) 
roughly agree with results by lChoyke fc Patrick! ( 119571 ) for colorless samples, but their values 
are lower than extrapolation of our difference results for the tan platelets or with the trends 
for the yellow and green slabs. 

A steep rise to the UV for our deeply colored gr een and yellow 6H SiC begins near 
21,000 cm~^ whereas spectra for the yellow 3C sample of iPhilipp &: TaftI (Il960l ) upturns near 
32,000 cm~\ and spectra for their colorless samples upturn at 35,000-38,000 cm^^ The 
previous absorbance data were corrected for surface reflections using n from microscopy. 
Although R values from microscopy are slightly larger than the direct measurements (Fig. [8]), 
small differences in R give a negligible logarithmic correction term and are not important 
in deterrnining A. Also, the microscopy data were collected from colored SiC samples by 
Thibaultl (119441 ). which have higher A and R than the colorless variety. The platelets are 
actually fairly dark, because color is apparent at 25 /im thicknesses, so it may not be valid to 
make a direct comparison to the previous data obtained from light colored samples (i.e. our 
darker samples contain more or different impurities). Given the steep rise for all samples, 
our linear extrapolation probably underestimated A and overestimated R above 31,000 cm~^ 
(Fig. [8]), and the high-order polynomial fits are probably the better representation (Fig. [8]). 
Nonetheless, the extrapolation in A is consisten t with that obta i ned fr om the Kramers-Kronig 
analysis of UV reflectance data for 6H SiC by lPhilipp fc TaftI (119601 ). 



Note also that n of 3C equals that of 6H SiC when light is polarized perpendicular to 
the c crystal axis; this is consistent with behavior at longer wavelengths. All evidence points 
to impurities, not structure, being the main cause of variation in the visible-UV region, as 
expected because electronic transitions are excited at such high frequencies. 

The difference method (§[3]) calculates a decrease in R between u = 23,000 and 29,000 cm~^, 
which is not in accord with classical dispersion analysis. This behavior is an artifact and 
is possibly due to the impurity content of the two platelets differing. It is more likely that 
the surface smoothness differed slightly. Variation in the amount of scattering is the most 
likely explanation for the dip in R. The steep rise in R, however, is attributed to impurities. 
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given the correlation in the upturn in A with the intensity of coloration. Reflectance will 
also upturn as R and A are correlated. We propose that the impure samples have more 
band structure in the UV than does pure, stoichiometric SiC, and reflectance increases at 
lower frequencies for deep colors. The linear extrapolation underestimates A and overesti- 
mates R. This behavior and the trend in n from microscopy of colored SiC suggests that 
impure, colored SiC samples have either an additional lower lying UV reflectivity band than 
the colorless variety, or perhaps that the entire band complex (Fig. [9]) is shifted to lower 
frequency. 

The similar raw absorbance obtained for both platelets and a much thicker moissan- 
ite sample (Fig. [7]) is consistent with absorption coefficient being essentially negligible in 
the near-IR to visible region. For this case, limiting values (see § |3]) are a better means of 
constraining optical functions than the difference method because errors result from sub- 
traction of measured values of the absorbance that are very close to one another. First, 
absorbance and reflection measurements of the 2.91mm moissanite crystal constrain A and 
R near 4000 cm~^ (Fig. [8]). Second, additional constraints are provided by the very thick 
green and yellow 6H slabs (Fig. [7j). The green 6H sample has a strong absorption feature near 
u = 16,000 cm~^ and the yellow 6H slab absorbs weakly in this spectral region. Absorption 
coefficients were calculated for the green and yellow 6H slabs by using R = 0.2 (Fig. 8, also 
see Part I), which is reasonable because R varies little from u = 9000 to 23,000 cm~^. Agree- 
ment is good with the difference calculation (Fig. [8l), especially below 16,000 cm~^ where the 
Bomem system operates well. 



4-3.2. Classical dispersion analysis of previous UV reflectivity data 



Reflection measurements of 6H by iPhilipp fc TaftI (Il960l ) are consistent with a grad- 
ual increase i n R from mois sanite measurements in the near JR. Kramers-Kronig analysis 
performed by iPhilippI (119581 ) gives A values that ar e consis t ent w ith our difference method 
but that are much larger than the measurements of iPhilippI (119581 ). Because of this discrep- 
ancy, we performed a classical dispersion analysis (see Part I) of iPhilipp fc TaftI (Il960l )'s 
reflectivity data. Using flve oscillators with positions, FWHM and oscillator strengths (/) 
shown in Fig. [9l we were able to reproduce their measurements as well as our values for R of 
moissanite in the near JR. We found that the best fits were obtained by setting = 1.56. 
Measurements of n are made in the visible region and do not offer an independent constraint. 
Although good agreement was obtained with Kramers-Kronig analysis at low frequency for 
n, the results for n are slightly lower than direct measurements of n by microscopy (Fig. [8]). 
Values for k do not agree. When overlapping peaks are present, as is the case here, Kramers- 
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Kronig analysis can fail, but the appearance of the optical functions does not suggest this 
is the case. Generally speaking, classical dispersion analysis is considered more trustworthy 
because it does not require extrapolation of R from zero to infinity. The classical dispersion 
analysis provides much higher values of A (Fig. [H]) than even the difference method. Such 
high values are not correct given the available absorption data, and point to small errors 
existing in the absolute measurements of R at low frequency, i.e., where the detector was 
least sensitive. 

As discussed in § 14.2.11 stacking disorder occurs for hexagonal varieties of SiC. This 
presents a potential problem when measuring R of 6H, because combined with small sample 
size it can lead to back reflections. Whichever of these problems exists is manifest at low 
frequency, where A is lowest. Imperfections in surface polish are likely for the platelets. This 
inflates R but not A in the difference method. Errors in R in the UV may also be instrumental 
effects for these very difficult, vacuum measurements. Absolute values for reflectivity also 
depend on the standard used for comparison (generally metallic silver). 



4-3.3. Construction of high frequency optical functions for SiC 



This section combines the available data to provide two sets of high frequency optical 
functions for SiC, one that represents stoichiometric, pure SiC, and another that represents 
impure SiC. Both sets are valid for the 3C structure and 6H with EJ^c. 

The starting point at 4000 cm~^ for pure and impure SiC is provided by our measure- 
ments of k and n for the 2.91 mm thick sample of moissanite. For impure SiC, electronic 
Table [7] is a linear interpolatation of k from the moissanite value at 4000 cm~^ to k for the 
platelets at 20,000 cm~^ obtained from applying the difference method to the low order poly- 
n omial and incorporates v alues obtained from A as shown in Fig. [8] to 30,000 cm~^. The fit 
of IShaffer &: Naunj (119691 ) constrains n from 10,000 to 30,000 cm~^ with linear interpolation 
below. Above 50,000 cm~^. Table [7] incorporates the classical dispersion analysis for n and 
k. Between 30,000 and 50,000 cm~^, linear interpolations were made. 

For pure SiC in electronic Tabled k betwe en i> = 36,000 to 56,000 cm~^ is constrained 
using the absorbance of 6H a-SiC measured by iPhilippI (Il958l ). whereas classical dispersion 
analysis is used above 70,000 cm~^, and interpolation is made between the segments. For 
n > 10,000 cm~^, classical dispersion analysis is used , and interpolat i on to moissanite data 
is made below 4000 cm~^. Classical dispersion fits to lPhilipp fc TaftI (Il960l ) UV reflectance 
measurements of 6H provides for extrapolation to infinity. For more accurate characterization 
of both impure and pure SiC, additional measurements of n or of i? from gemmy, colorless 
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moissanite are needed for both polarizations in the visible and UV. 



Discussion 



In carbon-rich dusty environments, SiC is one of only a few abundant dust species that 
exhibits an infrare d spectral feature from which the physical parameters of these regions can 
be diagnosed (see ISpeck et al.l l2009l . and references therein). Consequently, SiC has been 
of great interest to researchers seeking to understand the evolution of the dust shells and 
infrared features of carbon stars and their successor, C-rich pre-planetary and planetary 
nebulae. Furthermore, the lack of observable SiC spectral f eatures in the ISM has lead to 
many studies which place upper limits on its abundance (e.g.. lWhittet et al.lll99Cll : iMin et al. 



20071 : iKemper et al.ll2004j ). Moreover, the occurrence of presolar SiC grains provides the 
link between dust production in aging stars (AGB stars and supernovae) and subsequent 
formation of planetary systems, and thus SiC serves as a tool for dust studies in many 
environments. 

Astronomical studies of C-rich environments utilize two forms of laboratory data: (1) 
absorbance spectra, which are often converted to mass absorption coefficients for direct com- 
parison to observational spectra, under the assumption that the grain-sizes and shapes of the 
laboratory samples are similar to those dispersed in space, and (2) optical functions, which 
are use in various radiative transfer models and in analyzing the effect of grain morpholo- 
gies on spectral features, in which case single-reflection events are commonly assumed. In 
§ SHU we presented absorption coefficients for single- crystal SiC samples with both 3C and 
6H structural types and different amounts of impurities measured over a broader wavelength 
range. We focused on near-IR and visible frequencies because this region has been largely 
neglected. Merging our IR to UV measurements with previous far-UV data allowed us to 
provide continuous wavelength coverage for absorbance and optical functions of SiC, thereby 
permitting better constraints on its extinction for all astroph ysical environ ments. To pro- 
vide n and k over all wavelengths, we expanded upon work by iLipsoru (jl960l ) and developed 
a method to extract optical functions from transmission measurements which to the best 
of our knowledge has not been used heretofore in astronomy. The basis is that reflection 
processes are an intrinsic part of transmission/absorption experiments (Fig. [J). In Section 4 , 
we compared our new data with previous measurements at high (e.g., iPhilipp fc Taftlll960l ) 
and low (Part I) frequency, and found that values for optical functions are indeed high in the 
infrared. Through the comparison, we delineated the strong effects of polarization and im- 
purities. This section discusses limitations of the method to derive n and k and astronomical 
implications of the method and our results. 
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5.1. Caveats in using difference method and limiting values to obtain n, k 

The difference method can be used to calculate n and k for any material from absorption 
measurements of crystals with different thicknesses. However, accurate values are obtained 
only under certain conditions. Except for thickness, the samples used must be identical: the 
orientation and impurity content must be the same, and the surface polish must be uniform 
for all samples. Imperfections in the surface polish do not seem to affect A but alter R 
from intrinsic values; this affects the calculated n strongly and A; to a minor extent. In 
addition, the measured absorbance cannot be minuscule. For thin samples where A is near 
zero, limiting values provide better constraints. Measuring absorbance of additional, larger 
crystals and obtaining independent measurements of the n will constrain n and k at all 
frequencies above the main peak. Using large, single- crystals is the most accurate means of 
determining k for wavelength regions in which k is low. The difference method is particularly 
useful in the visible region, where absorbance is high but k is low (see Eq. [51 ) 

We found that the best approach is to collect data from single- crystals of diverse thick- 
nesses and to use both limiting values and the difference method. For SiC, our analysis 
was complicated by the variation in color (impurities) for our samples and those previously 
studied. Color differences clearly alter absorption coefficients in the visible-UV region. We 
provide end-member values for n and k, based on presence and absence of color. The differ- 
ence methods and limiting values cannot be applied to measurements of powder dispersions 
(see §[E2l). 

5.2. Why A;(A) derived from thin film, powder dispersion, reflectivity 

measurements differ 

Although k is a material property and should in principle be independent of the method 
used to determine it, in practice, values of k for any given substance are influenced by the 
experimental technique. The underlying problem is that laboratory studies do not directly 
measure absorption, but rather determine transmission which contains a reflection compo- 
nent, and, moreover, at any given wavelength, R and A are correlated. Near the frequencies 
where light is strongly absorbed, light is also strongly reflected. For SiC in particular, R 
approaches 97% near A of 10 fim. Thus, at frequencies near the SiC fundamental mode, 
large particles are effectively opaque and the 3% of light that is transmitted is attenuated 
sufficiently that the small amount which emerges from the sample is within the noise level 
of detection. For SiC, particles with sizes near 1 fim are effectively opaque. This size is con- 
sistent with grains in space and in powder dispersions used in laboratory studies. The effect 
creates non-linear dependence of transmission measurements, used to determine absorbance. 
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on grain size for strong IR modes. For small-grain dipsersions (e.g., in a KBr pellet), the 
grains can clump, increasing the effective grain size to larger than the average created during 
grinding. This skews the grain-size distribution and makes a greater portion of the particles 
in the pellet effectively opaque, which is problematic as optically thin conditions are assumed 
in analyzing laboratory spectra. 

Second, in addition to the processes of reflection and absorption, a non-negligible amount 
light is transmitted unattenuated through parts of the pellet matrix; powder dispersions 
transmit some unattenuated hght at all frequencies because the particles do not completely 
cover the area. Such light leakage makes the fundamental peak in dispersion spectra of 
SiC both more rounded and lo wer in measured A or de rived k values at peak center when 
compared to thin film spectra (ISpeck et al.l Il999l . 120051) or in ideal absorpt i on co efficients 
(Part I). For detailed disc ussion of light leakage (see iHofmeister et al.ll200ll . l2003l ). SiC in 
particular is discussed by ( Hofmeister et al. 200ol ). 



The two effects combined mean that A and thus k (Eq. [5]) obtained from powder dis- 
persion methods is lower than the bulk grain-size independent property near the main peak 
obtained from reflectivity data. Measurements of k from thin films are a little less than from 
reflectivity data, particularly at the main peak, but larger than k from powder dispersions, 
because light leakage exists in the film, but is much less than in the dispersion. In essence, 
it is the presence of reflectivity in all types of absorption measurements that causes the 
variation in the derived values of k from intrinsic, reflectivity measurements. 

These differences between measurement techniques preclude the use of powder disper- 
sions with the difference method because (1) the thickness of the material in a pellet sample 
is not well constrained; and (2) multiple reflections may occur from the particulates in the 
powder. 

Consequ ently, thin film absorbance sp ectra offer an improvement over powder dispersion 
spectra (e.g.. IHofmeister et al.ll200lU2003l ). but because film thickness is difficult to measure 
accurately, quantities such as k, band strengths, and extinction coefficients are accompanied 
with experimental uncertainties. In this paper, we constrained thickness of films previously 
used to collect mid-IR data through comparison with near-IR spectra of single crystals. The 
resulting high values of k for the main Si-C stretching mode support the results of Part I 

The presence of reflections and light leakage in transmission measurements makes deriva- 
tion of mass absorption coefficients from powder dispersion data highly uncertain. This 
problem is exacerbated by the potential effect of grain shape. Many astronomical studies 
use mass absorption coefficients obtained from powder dispersion data, /tabs, to estimate 
the relative abundances and total masses of dust species in a given circumstellar shell (e.g.. 
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Albrecht fc Chinil I2OO0I : Istevens fc Geai] I2OO0I : iKemper et all I200I booj ). In light of the 



hmitations on using powder dispersions, the vahdity of these derived abundances is suspect. 

Using the new optical functions presented here, it is possible to calculate /tabs for arbi- 
trary grain sizes and grain shapes. Furthermore, the effect (or lack thereof) of impurities 
can be assessed. Thus the errors on previous derivations of abundances and dust masses can 
be estimated and true limits of these properties of dusty regions can be achieved. 



5.3. Importance of reflectivity in modeling extinction in the near-IR to visible 

regions 

In radiative transfer models, light that is extinguished by dust contains absorbing and 
scattering components. 



Qext Qahs ~l~ Qsca 

This section shows how the Q factors are related to material properties for the spectral 
regions focused on in the present paper and ties our results to radiative transfer modeling. 

The exact method by which Q-values are calculated varies, and the simplest case (spher- 
ical grains) uses Mie th eory to compute th ese efficiency factors from the complex index of 
refraction m = n — ik (Ivan de Hulstlll98ll ). It is unlikely that astrophysical grains are all 
spher ical, and recent studies have investigated the effects of grain shapes on spectral feature s 
(e.g. , ISohren fc HuffmanlEosal : iMutschke etallll999l : iMin et allBoOsI : kndersen et al.ll2006[ ). 



However, the underlying relationship between the scattering and absorption efficiencies {Q 
values) and laboratory measurements of i?, A from which Qs are derived will apply to all 
grain shapes. Here we use Mie theory as an illustration. 



From Mie theory: 



Qahs — 4x55 



m 



(10) 



^ 3 \m^ + 2 



where x = 27r/i/A and h is grain size. It is assumed that x <^1. 
Using the definition X = l/u and Eq. [5] we obtain: 

Qahs = A 



(11) 



12nh 



F + 2)2 + 4n2F 



(12) 
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From Eqn. 5 and Figures 8 and 9, A; is near zero for the transmitting near-IR to visible 
region. In this region n is small but finite and not very frequency dependent. Consequently, 
the denominator in Eq. 12 is near the constant 4, and the numerator is proportional to the 
product An. Through this relationship between Qabs and A, the low values of absorptivity 
control values of Qabs making it near zero in the near-IR to visible region. Thus, extinction 
is due to scattering, not absorption, over this particular spectral region. 

To further probe the connection of radiative transfer with material properties, we sim- 
phfy the expression for Qsca.'- 

Qsca=l{2^huY{^^' (13) 

This is proportional to R from eq. [3 

g.,, = ^-{2^huYR {^^^^ (14) 

Again, in the near-IR to visible spectral region, n is small but finite and weakly dependent 
on frequency. To first order, the numerator and denominator involving n cancel, and Qsca in 
the near-IR to visible region is controlled by its strong dependence on frequency and by the 
specific values of R. Our analysis shows that refiection is important in radiative transfer, not 
only in spectral regions where absorption dominates, but even more so in spectral regions 
characterized by weak to negligible absorptions. Constraining R in the visible through 
laboratory measurements is important. 



5.4. Potential uses of new data on n and k for SiC over the IR to UV 



In Part I we provided quantitative mid-IR laboratory measurements of SiC, which al- 
ready improved the data available for astronomers in a variety of dusty fields, most notably by 
including /3— SiC, rather than the a— polytype, and by having truly grain-size independent 
data. However, the new data presented here augments Part I by extending the wavelength 
coverage blueward. This is important because SiC grains are found in environments where 
the radiation field is strong in visible or UV regions. Furthermore, our data show that while 
the main ll.Syum peak may not be diagnostic of grain impurities, the UV absorption is. 

In addition, it is n ow possible to use the SiC optical funct ions to analyze grain-size and 
grain- shape effects (c.f. iBagnulo et al.lll995l : ISpeck et al.ll2005l ). It has been suggested tha t 
there is an evolution in SiC grain sizes with the evolution of the star fjSpeck et al.ll2005l ). 
Studies of the effect of grain size evolution have not been possible not only because of a 
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lack of grain size independent data, but also because as the host stars evolve, the radiation 
field contains more high energy photons and the optical data needed for modeling has been 
lacking. Furthermore, the SiC feature prevalent in carbon star spectra becomes much less 
common as these stars progress to the end of the AGB phase and into the post-AGB and 
PNe phases. 

Whereas SiC is commonly seen in C-s tar spectra, its occurrence in Galac t ic (Milky Way) 



plane tary nebulae (PNe) is less common (lBarlowlll983l : iRinehart et al.ll2002l : ICasassus et al. 



2OOII ). even though C-stars evolve into C-rich PNe. However, evolution of these grains from 
the C-star to PNe phase is not understood, and other spe ctral features (e.g., ~ 11 — 12/im 
plateau due to PAH emission) may mask the SiC band (see lSpeck &: Hofmeisterl 120041 ). Even 
though SiC is apparently rare in Galactic PNe - there are a handful of examples. We need 
to include SiC optical functions into the models to see how much can be hidden. 

This trend in the occurrence of the SiC feature from C-stars to pre-planetary nebulae 
and PNe is not echoed in recent studies of the LMC. The SiC feature remains observable, 
and thus including SiC into models of Magellanic cloud objects is essential. Knowledge of 
optical properties of SiC is essential to understanding dust formation in lo wer metallicity 



environments such as the Mage l lanic Clouds (MCs ) and local group g alaxies (ILagadec et al. 



20071 . 12OO8I : IZiilstra et all I2OO6I : ISpeck et all I2OO6I : ISloan et al.ll2008f ). Furthermore, recent 
observations of PNe in the Magellanic C louds show that SiC features are much more commo n 
in these low metallicity environments (jStanghellini et al.l 120071 : iBernard-Salas et al.ll2008l ). 
In addition, the recent discovery of extreme (optically-obscured) C-stars in the LMC shows 
that SiC absorption features also occur more often in these low metallicity environments 
(IGruendl et al.ll2008l ). Understanding why SiC features are more common in these highly 
evolved carbon-rich environments than their Galactic counterparts requires reliable, grain- 
size-independent optical functions for SiC. Moreover, modeling of PNe requires wavelength 
coverage into the ultraviolet (UV) region, as this is where the radiation from the exciting 
central star peaks. 

The fact that we can distinguish between samples of different purities based on the 
UV extinction has potentially far reaching consequen ces. In particular, t he inclusion of 
SiC into dusty photo-i onization codes like MOCASSIN jErcolano et al.lboosh and CLOUDY 
(jVan Hoof et al.ll2004j ) is desirable. It will be very interesting to see the impact of varying SiC 
purity on the overall outputs of these codes. Furthermore, the ability to use UV observations 
to determi ne the purity of the SiC g rains provides a valuable test of the hypotheses put 
forward by lSpeck fc Hofmeisten (120041 ) on the "21/im" feature, for example. 



The changes in the UV absorption properties could potentially have huge impacts on 
hydrodynamic models of PNe shaping where the details of the interaction of light and dust 
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grains is important. 



6. Conclusions 



We have discussed in detail the relationships between reflectance and transmittance 
measurements and how material properties such as optical functions can be extracted from 
such data. We emphasized that absorbance and absorption coefficients are derived, not 
measured, quantities, and that reflection plays a crucial role in both laboratory measurements 
and radiative transfer models. 

Although absorbance can be compared directly to observational spectra in the IR, this 
does not take the effect of grain shape into account. The extracted n and k values can be 
used for more detailed analyses of astrophysical dust. 

Reliable, grain-size/shape independent optical functions covering a broad wavelength 
range for a range of minerals is a growing need in the astronomy community. To address 
these needs. 



1. we have devised a new "difference" method, in order to obtain n and k from single- 
crystal absorption data. For spectral regions where both reflectivity and absorptivity 
are low but non- negligible, transmission data from two sections with different, known 
thicknesses reflect the same amount of light at the front and back surfaces, while 
simultaneously absorbing different amounts of light. The transmission data thus yield 
two equations with two unknowns, which can be solved uniquely for n and k. If used 
carefully, the difference method and use of limiting values can provide n and k over 
a broad wavelength range for astronomically relevant minerals using widely available 
transmission techniques. 

2. We presented new SiC absorbance data which both extends the coverage of existing 
data from IR to UV wavelengths and checks optical functions deriv ed from reflectance 



measurements. To improve upon our previous thin-fllm IR data (ISpeck et al.l Il999 



20051 ) . we merged it with new single-cr ystal spectra at slightly h igher frequencies 



thereby constraining fllm thickness (after iHofmeister fc Boweyl l2006l ) and allowing ac- 
curate extraction of k. This also r nakes direct co r npari son with absorption coefficients 
extracted from reflectivity data in iPitman et al.l (120081 ) possible. 



3. We implemented the difference method to provide IR to UV optical constants for SiC 
[v > 1000 cm^^) for 611 SiC with EA^c. These data are also valid for 3C SiC, the (3 



- 25 - 



structural type more dominant in astrophysical environments but are affected by impu- 
rity content at high frequency. Infrared peak positions for ii^||c are shifted from those 
of 3C, particularly for the overtones, permitting the 6H species to be distinguished. We 
have thereby extended our previous reflectivity determinations of n and k for SiC into 
the UV spectral region. We also fitted available UV reflectivity to Lorentz oscillators 
as another check and to extend the frequency range. 

4. We provided detailed measurements on SiC overtone-combination bands, which are 
affected by structure and orientation. We have shown that the impure (colored) and 
pure varieties of commercially prepared SiC produce identical k values just below the 
visible. We have also confirmed that the 3C spectrum is n e arly i dentical to that of 6H 



when oriented with El.c as determined in iPitman et al.l (120081 ) for the fundamental 



frequency modes. Available data indicate that the equivalence holds for all spectral 
regions, if sample purity is the same. 



Our new approach - the difference method - for extracting optical functions from trans- 
mission spectra can be applied to any material, not just SiC. Because transmission is easily 
and commonly measured, our difference method may be used to greatly expand the database 
of astronomically relevant near-IR to UV optical functions and aid implementation of radia- 
tive transfer models for many environments, not just C-stars. 
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Fig. 1. — Schematic of light loss in a dielectric material showing the most intense ray 
paths at near-normal incidence (the angle of incidence is exaggerated here to provide a clear 
representation of the processes). The incident light with intensity Iq is partially reflected and 
transmitted at the front surface. The transmitted beam is partially absorbed, determined 
by thickness d, reflectance R, and absorptivity = fl = I^hsl h- At the back surface, the beam 
attenuated to /o(l--R)(l-f2) is again partially reflected and transmitted. For the case of a very 
weakly absorbing and thin sample, and if the angle of incidence allows the beam to reach 
the detector, then Igum = 2RIo approximates the amount of light measured in a reflectance 
experiment. In contrast, R is typically 20% for SiC in spectral regio ns other than the naid-IR , 
and should not be neglected in analyzing transmission data. After iHofmeister et al.l (120031 ). 



- 31 - 



Wavelength, |Lim 
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Fig. 2. — Raw absorption spectra of /3-SiC. Values are not corrected for any baseline or for 
reflectivity losses. Thin solid line = 3C wafer with thickness = 308 fim. Thick dotted line = 
same sample at cryogenic temperatures. Thick solid line = same wafer, thinned to 25 fim. 
Dashed line = thin film of nano crystalline /3-SiC from Speck et al. (2005): only this sample 
pertains to the left side axis. Interference fringes exist in this spectrum. Vertical lines denote 
weak peaks with positions given in wavenumbers. 
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Fig. 3. — Baseline-corrected and merged absorption spectra of /9-SiC. Thin dotted line = 
merge to thin-film spectrum, scaled to match at high frequency. The merge suggests an 
original film thickness of 0.74 /zm. Thin solid line = merge to thin-film scaled to match 
at low frequency. The merge suggests an original film thickness of 0.34 /zm. Thick solid 
line = ideal absorption coefficient calculated from classical dispersion analysis of measured 
reflectivity. 
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Fig. 4. — Absorption spectra of 6H SiC, E-Lc, from A l fa/Ae sar. Except for the thin-film 
made from crushed platelets presented in ISpeck et al.l (120051 ) . all spectra are from single- 
crystals. Thin line and gray line = very thin platelets with nominal thicknesses of 5 /xm. 
Heavy line = platelet with a nominal thickness of 15 /im which was scaled to match spectra 
of the 5 fim platelets. The scaling used suggests that the sample actually was originally 
22 fim thick. Dotted line = thin film. From scaling (light dotted line), the thickness of the 
film is 0.38 /im. Because interference fringes exist in most spectra, intrinsic peak positions 
are indicated by vertical lines with labels in cm~^. 
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Fig. 5. — Polarized absorption spectra of moissanite (6H-SiC). Broken lines = E±c. Solid 
lines = E\\c. For the thinnest sample (light lines), the measured thickness was 85 iim. 
The moderate thickness sample (heavy lines) was rougly 250 fim. thick. Thickness is well- 
controlled for the largest sample (medium lines on the right-hand side, corresponding to the 
right-hand scale.) 
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Fig. 6. — Merged polarized absorption spectra of moissanite at high frequency with the 
platelets at low frequency for E1.C. Thick dotted line = ideal absorption coefficient calcu- 
lated from classical dispersion analysis of measured reflectivity. Thick solid line = merged 
spectrum, scaled assuming platelet thickness of 5 /xm and film thickness of 0.18/im. Dashed 
line = merged spectrum, scaled assuming largest moissanite slab was 2.91 mm thick, i.e., 
film was 0.54 /xm thick. Fine solid line and right axis = expanded view of merge of moissanite 

— * 

at high frequency with classical dispersion analysis for E^c. Fine dotted line and right axis 

— * 

= merged spectrum for El.c. Overtones are scaled based on the measured thickness of 2.91 
mm. 
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Fig. 7. — Unpolarized near-IR to UV spectra from 6H platelets and colored slabs. Fine black 
curve = thicker platelet. Gray = thinner platelet. Heavy black curve = small moissanite 
gemstone. Dashed line = yellow 6H slab. Dotted line = green 6H slab. Thicknesses as 
labeled. 
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Fig. 8. — Properties derived from measurements in Fig. [7] and comparison to previous data, 
(a) A and R. Dotted line = true absorption coefficient obtained by difference. Heavy double 
arrow denotes region where R seems underestimated, see text. Light gray line = A for green 
6H SiC calculated by using constant R = 0.2. Dashed line = A for yellow 6H calculated 
by constant R of 0.2. Heavy dotted line near 4000 cm~^ = A from near-IR measurements 



Philippl fll958h. Circle s 



of the thickest mo issanite s lab. D iamonds = A measured for 3C 
= A measured by IPhilippI (1l958l ) . Open plus = A measured by 

Dot-dashed = A fr om classical dispersion analysis of R from by lPhilipp fc TaftI (jl960l ). Small 



Chovke fc Patrick! fll957l ). 



squares = A from lPhilipp fc TaftI (Il960l ) using Kramers- Kronig analysis of R. Solid line = 
reflectivity obtained by difference. Heavy solid line = reflectivity measured for moissanite 
(Part I). Heavy grey l ine = R from classical dispersion analysis of reflectivity data (squares) 
from iPhilipp fc TaftI (Il960l ) which is indistinguishible from their Kramers-Kronig analysis, 
(b) Optical functions, mostly for 6H. Dot-dashed line = k hj difference using a high-order 
polynomical fit over all frequencies measured. Heavy dotted line = k from near-IR mea- 
surements of moissanite. Light solid line = n hj difference using a high-order polynomical 
fit. Horizontal double arrow = region where polynomial fits were not perfect. Heavy line 
= n from near-IR refle ctivity data on moissanite of Part I. Grey line = n from microscopy 
( IShaffer fc Naumlll969l ). linearly extrapolated from the measured spectral region, denoted 
by arrows. Tr iangles filled = k and open = n from Kramers-Kronig analysis of reflectivity 
performed b v iPhilipp Sz TaftI (Il960l) . Dots = n from our classical dispersion flts to mea- 
surements of iPhilipp fc TaftI (Il960l ). Dashes = k from our classical dispersion analysis of 
Philipp fc TaftI (Il960l ) data. Circles = k from absorbance measurements of colorless 6H by 
PhihppI (ligsi l 
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Fig. 9. ^ — Classical dispersio n analysis (CDA) of UV reflectivity data from = El.c of 6H 
SiC by iPhilipp fc TaftI (119601 ) Jagged line at low frequency is UV reflectivity of moissanite. 
Squares = measured values of R. Short dashes = CDA results for R. Dotted curve = CDA 
results for n. Long dashes = CDA results for k, based on the peak parameters listed in the 
inserted table, v is peak position; FWHM is full width at half maximum, both in cm ~^; / 
is oscillator strength. Solid lines = Kramers- Kronig analysis by iPhilipp fc TaftI (Il960f) fthin 



for 77- , thick for k). Circles and dots are polarized measuren ients of n for 6 H SiC by iThibault 
(119441 ) Thick black line = k for the colorless 6H sample of jPhilippI (119581 ). 
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Table 1: Experimental samples descriptions 



Polytypc Mineral Name 



Grain Size 



Comments 



a-SiC synthetic moissanite 

a-SiC 6H-SiC (amber) 

a-SiC 6H-SiC (blue) 

a-SiC 6H-SiC (green) 

a-SiC 6H-SiC (yellow) 

/3-SiC 3C-SiC wafer (fee cubic) 

/3-SiC 3C-SiC gray (fee cubic) 



diam. = 6.5 mm 
2 yum powder; surf, area 
= 9-11 m^g-i 
several mm per crystal 
diam. ca. 4 mm 
diam. ca. 8 mm 
diam= 5 /im 
diam. ~ 2.5-25 /xm 



round brilliant cut gem 
hexagonal plates; purity: 99.8% met- 
als basis 

Intergrown crystals 
layered hexagonal crystals 
single crystal, some zoning 
CVD wafer purity: > 99.9995% 
equant (spherical) chips 
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Table 2: Frequencies of overtone-combination modes in SiC, in cm 

3C 6H 

Mode Type 





Calc. 


Meaa 


Calc.'^ 


Meas ^ 


o iA = i U 


l\) i 




7n7 
l\) i 


707 ^ 1 




OOi 


OOi 


oyi 




i u-|-roiaea 


n.a. 


1 ni 
iUiZ 


iUzZ, iUo/ 


1 nnn 
iUUU 




iUD4 


iUMD 


iU04,iUD4 


1 n7Q 1 


T A _1_0T A 


11/17 
ii4 ( 


1 1 Ot^V. 
iiZOD 


1 1 c;7 

iiO ( 


1 1 QAK 


9T A 


iZZO 


1 on7 

iZU < 


iZZO 


1 1 O/l 

iiy4 


Li(J-t- iA 


1 OQO 




1 OQO 

izoy 




1 U+zroiu 


n.a. 


iZo4 


1 070 1 
iZ ( Z _L 


1 OC*^ 1 
iZoO _L 


± (j-\-z lA 


1 QOO 

iozy 


1 Q1 1 o 
ioiiS 


ic)iy, iozy 


ioUoS 


7 








1 QQ7 1 
ioo / _L 


T^O_lT a 
1 U+1jA 


1 /1 1 Q 


i4UZS 


1 /1 1 Q 1 /10Q 
i4io, i4Zo 


1 /IH/I 1/11 Oa 
i4U4, i4iZ 


OT A 1 T^A 
ZijA+ iA 


1 /lO^i 

i4yD 




1 ^1 
iOiD 




T O-LTA 
Li(J-\- i A 


iOUO 


iOZOS 


iOUU, iOUO 


iOZU, iOZD 


7 




iOoOD 




i04o 


T 0_1_T A 


1 c;qq 
iooo 




it)oo,iooo 


1 err r -1 Cfjn 
iOOO,iODU 


Z i <J 


ioyo 


iOZOS 


io/0, ioyo 


1 A1 S 1 fJOHo 
iOio, iOZUS 


7 








iD4o 


TO+TA+LA 


1678 




1678, 1688 


1683, 1700 


2TA+2LA 


1762 


1716 


1782 


1713,1720b 


LO+TA+LA 


1854 




1860, 1864 




2T0+TA 


1861 


1920b 


1841, 1861 


1876, 1890b,s 


? 








1897 ||s 


T0+2LA 


2027 




2037, 2047 


1980, 1990 


2TO+2TA 


2128 


2079 


2108, 2128 


2076, 2080s 


2T0+LA 


2209 


2182b 


2210, 2219 


2179b 


T0+L0+2TA 


2302 


2297s 


2287, 2297 


2286, 2297s 


3T0 


2392 


2396 


2364, 2392 


2388, 2394s 


2T0+LA+TA 


2175 


2187 


2-1G5. 2185 


2177. 2-180 



Notes: b = broad, s = strong, n.a. = not applicable. "Calc" values based on fundamental 
frequencies (Part I, Nakashima & Harima 2000), summarized as follows in u (cm~^) (A in 
/xm = 10,000/z/): TA = 266 in 3C, also at 234.8 and 244 in 6H; LA = 614.4 in 3C, or 624.4 
in 6H; TO = 797.7 in 3C, 797.5 for E^c , 787.8 for E\\c- LO = 973 in 3C, 975 for i±c , 970 
for E\\c. 6H also has modes at 883.7 and 888.5 for E\\c. 

^ For 6H calculated, if two entries are given, the lower frequency value pertains to E\\c. Modes 
marked || or _L are found only when E is oriented parallel or perpendicular to the c-axis. 
Unmarked entries indicate that the same frequency exists in both polarizations. 

^ For the measured modes of 6H, modes in both polarizations are unmarked. If two entries are 
given, the lower frequency value pertains to E\\c; for these split modes, entries having lower 
values in £^||c confirm involvement of the TO or LO components. 
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Table 3: Merged single-crystal al)soil)aiice spectnuri of (3C) ,5-SiC* 



Wavcnumber 


Absorption Coefficient 




(;um~^) 


401.12 


-0.0054155 


402.09 


-0.12127 


403.05 


-0.071841 


404.02 


-0.042450 


404.98 


0.70815 


405.94 


0.035504 


406.91 


-0.057883 


407.87 


0.060957 


408.84 


0.027782 


409.80 


-0.002GC73 



*Table continued as electronic-only file. 
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Table 4: Merged absorbance spectrum of (a) 6H-SiC, polarization (ordinary ray)* 



Wavenumber 


Absorption Coefficient 


(cm~"^) 


(/im~^) 


3.0000 


0.00020004 


4.0000 


0.00020007 


5.0000 


0.00020012 


6.0000 


0.00020017 


7.0000 


0.00020023 


8.0000 


0.00020030 


9.0000 


0.00020038 


10.000 


0.00020047 


11.000 


0.00020056 


12.000 


0.00020067 



*Table continued as electronic-only file. 



-43- 



Table 5: Ideal absorbance spectrum of (a) 6H-SiC, calculated from classical dispersion anal- 
yses, £^_Lc polarization (ordinary ray)* 



Wavenumber 


Absorption Coefficient 


(cm~^) 




3.0000 


0.00020004 


4.0000 


0.00020007 


5.0000 


0.00020012 


6.0000 


0.00020017 


7.0000 


0.00020023 


8.0000 


0.00020030 


9.0000 


0.00020038 


10.000 


0.00020047 


11.000 


0.00020056 


12.000 


0.00020067 



*Table continued as electronic-only file. 
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Table 6: Ideal absorbance spectrum of (a) 6H-SiC, ^||c polarization (extraordinary ray)* 



Wavenumber 


Absorption Coefficient 


(cm~"^) 


(/im~^) 


2.0000 


3.0000x10"^ 


3.0000 


7.0000x10"^ 


4.0000 


1.2000x10"^ 


5.0000 


1.8000x10-^ 


6.0000 


2.7000x10-^ 


7.0000 


3.6000x10-7 


8.0000 


4.7000x10-'^ 


9.0000 


6.0000x10-7 


10.000 


7.4000x10-7 


11.000 


9.0000x10-7 



*Table continued as electronic-only file. 
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Table 7: Complex refractive indices for impure SiC, applicable for either 3C (/?) or 6H (a) 
£^_Lc polarization (ordinary ray)* 



Wavenumber 
(cm~^) 


n 


k 




4000.0 


2.5574 


6.0687x10" 


-5 


4001.0 


2.5574 


6.0676x10" 


-5 


4051.0 


2.5575 


6.0984x10" 


-5 


4101.0 


2.5576 


6.1292x10" 


-5 


4151.0 


2.5578 


6.1599x10" 


-5 


4201.0 


2.5579 


6.1907x10" 


-5 


4251.0 


2.5580 


6.2215x10" 


-5 


4301.0 


2.5581 


6.2523x10" 


-5 


4351.0 


2.5583 


6.2830x10" 


-5 


4401.0 


2.5584 


6.3138x10" 


-5 



*Table continued as electronic-only file. 
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Table 8: Complex refractive indices for pure SiC, applicable for either 3C (/?) or 6H (a) EJ^c 
polarization (ordinary ray)* 



Wavenumber 
(cm~^) 


n 


k 




4001.0 


2.5501 


6.0687x10- 


-5 


4051.0 


2.5508 


6.0591x10- 


-5 


4101.0 


2.5515 


6.0495x10" 


-5 


4151.0 


2.5522 


6.0400x10" 


-5 


4201.0 


2.5529 


6.0304x10" 


-5 


4251.0 


2.5535 


6.0209x10- 


-5 


4301.0 


2.5542 


6.0113x10- 


-5 


4351.0 


2.5548 


6.0017x10- 


-5 


4401.0 


2.5553 


5.9922x10- 


-5 


4451.0 


2.5559 


5.9826x10- 


-5 



* Table continued as electronic-only file. 



